arXiv:1504.03776vl [quant-ph] 15 Apr 2015 


On the effects of self- and cross-phase modulation on photon purity for four-wave 

mixing photon-pair sources 

Bryn Bell/’Alex McMillan/ Will McCutcheon/ and John Rarity^ 

^Department of Eleetrieal and Eleetronie Engineering, University of Bristol, Bristol BS8 lUB, UK 
^Centre for Ultrahigh bandwidth Deviees for Optieal Systems (CUDOS), 

Institute of Photonies and Optieal Seienee (IPOS), 

Sehool of Physies, University of Sydney, NSW 2006, Australia 
(Dated: April 16, 2015) 

We consider the effect of self-phase modulation and cross-phase modulation on the joint spectral 
amplitude of photon pairs generated by spontaneous four-wave mixing. In particular, the purity of a 
heralded photon from a pair is considered, in the context of schemes that aim to maximise the purity 
and minimise correlation in the joint spectral amplitude using birefringent phase-matching and short 
pump pulses. We find that non-linear phase modulation effects will be detrimental, and will limit 
the quantum interference visibility that can be achieved at a given generation rate. An approximate 
expression for the joint spectral amplitude with phase modulation is found by considering the group 
velocity walk-off between each photon and the pump, but neglecting the group-velocity dispersion at 
each wavelength. The group-velocity dispersion can also be included with a numerical calculation, 
and it is shown that it only has a small effect on the purity for the realistic parameters considered. 
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INTRODUCTION 

Single photon sources are a vital component of de¬ 
veloping quantum technologies such as quantum cryp¬ 
tography [T], linear optical quantum computing [2], and 
quantum metrology [3], and improved sources need to 
be developed to enable these applications. In addition 
to requiring high efficiency, on-demand single photons, 
many applications require that the photons be generated 
in a single-mode, with well defined spatial characteristics 
and a fourier transform limited spectral-temporal shape. 
This allows two separate photons to be indistinguishable 
and to undergo high quality quantum interference [4] - 
this, in turn, makes possible fundamental operations such 
as teleportation of the photon [5], and two photon logic 
gates [6]. 

Photon pairs generated in a nonlinear medium by 
spontaneous parametric downconversion (SPDC) or four- 
wave mixing (FWM) are often used as a source of 
single photons, with one of the photons detected to 
give a heralding signal, indicating the presence of the 
other [3 i]. Although this method is inherently non- 
deterministic, through the multiplexing of many such 
sources and the use of active switching, it is in princi¬ 
ple possible to construct a source arbitrarily close to a 
deterministic source [3 HO]. However, the photons of a 
pair are generally correlated in frequency or time, which 
means that the single photons will arrive in a statistical 
mix of multiple modes. Narrow spectral filtering of the 
single photons can force them into a single mode, but 
at a cost to the overall transmission and heralding effi¬ 
ciency, which reduces the usefulness of the source. Pos¬ 
sible solutions to this problem have been demonstrated 
based on consideration of the joint spectral amplitude 


(JSA) of a pair: with careful choice of wavelengths, or 
by engineering the dispersion properties of the nonlin¬ 
ear medium, the degree of correlation can be minimised, 
allowing quantum interference to take place without nar¬ 
row filtering m- For SPDC in bulk crystals interference 
visibilities as high as 94.5% have been observed without 
filtering m- For FWM in birefringent optical fibres un¬ 
filtered visibilities have tended to fall short of theoretical 
estimates, often in the range 70 — 80%, which is far from 
sufficient for scalable use in communications or comput¬ 
ing [T3HT6]. It has been suggested that inhomogeneity 
along the length of fibres due to fabrication imperfec¬ 
tions is responsible for the short-fall in visibility, and it 
has been shown theoretically that a large inhomogenity 
can reduce the visibility HU- However, it is also expected 
to create a broadening and modulation in the spectra of 
the individual photons, which should be easily detected, 
and in some situations a small amount of inhomogeneity 
could actually improve the interference visibility m- 

Here, we consider the effects of parasitic nonlinear pro¬ 
cesses on the JSA and the interference visibility, namely 
self-phase modulation (SPM) and cross-phase modula¬ 
tion (XPM), which are not present in a nonlinear 
medium such as the crystals used for SPDC, but are 
potentially significant in a nonlinearity such as fi¬ 
bre m- Previous calculations of the JSA have tended to 
account for SPM and XPM in a simplistic form which is 
only exact in a continuous wave (CW) regime, whereas 
the relevant schemes to make the JSA uncorrelated rely 
on the use of short pulses. We show that these effects 
can cause a reduction in interference visibility, especially 
when a source is operated at a high pair generation rate, 
beginning from an analytic explanation then progressing 
to a numerical model taking full account of SPM, XPM, 
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and dispersion. 


JOINT SPECTRAL AMPLITUDES 

In pair production through FWM, a bright pump laser 
is used to power the process. As the pump pulse propa¬ 
gates through a medium, two pump photons may be 
spontaneously annihilated, with a correlated signal-idler 
photon pair created. The frequencies of the signal and 
idler are constrained by the conservation of energy and 
momentum: 


Auj = 2ujp — ujg — oJi = 0 .X 

Af3 = 2f3p - - A - 27 P = 0. ^ ^ 

ujp^s,i refer to the frequency of pump, signal, idler, and 
Pp,s,i fo the wave vectors. The 27 P term arises from 
SPM and XPM, as the intense pump in the nonlinear 
medium will slightly modify the wave vectors, with 7 an 
effective nonlinear coefficient, and P the peak power of 
the pump. 

These conditions are not exact and the photons will 
have some bandwidth centered on an exact solution. The 
finite bandwidth of the pump will introduce some uncer¬ 
tainty to the Auj condition, and for a fibre of finite length, 
there is some uncertainty in the phase matching which 
permits small values of Ap. The JSA can be simply 
expressed as the product of an energy matching and a 
phase-matching function m- 

JSA{u}s,u}i) = F X G ( 2 ) 


with 

F = jj (kjpidiAjp2E{iAjpi)E{u)p2)5{u)pi + ujp2 - Wg - w*) 

(3) 

G = (4) 

Here, E{uj) is the spectral amplitude of the pump, and 
the two pump photon frequencies ujpi and ujp 2 are inte¬ 
grated over. It can usually be assumed that the two are 
approximately equal, then F is just the convolution of 
E{uj) with itself, and the delta function, which results 
from energy conservation, fixes ujp = {oUg -\-uJi)/2. For 
instance if the pump amplitude is a Gaussian, E{ujp) = 

Eqc 2^2 ^ then 


(<^s+<^1 — 2wpQ)^ 

F = EqC 16ct2 


(5) 


Another useful simplification is to consider the differ¬ 
ent group velocities at the pump, signal, and idler fre¬ 
quencies but ignore higher-order dispersion terms. Then, 



FIG. 1. Joint spectral amplitudes calculated simply from the 
pump, signal, and idler group velocities, (a)General example 
with l3ip > Pu > Pis results in a highly correlated signal 
and idler, indicated by the diagonal nature of the main peak 
(b) Pip = Pii results in the central peak becoming vertical 
and uncorrelated (c) when Pip — Pis ~ Pu — Pip^ the phase¬ 
matching condition lies at +45° while the energy matching 
lies at —45%. Tuning the pump bandwidth can make the 
central peak circular and uncorrelated. All JSAs are plotted 
as absolute values, ignoring the complex phase. 


for small departures from an exactly phase-matched solu¬ 
tion Aujg and AuJi, the phase mismatch can be expressed 
as: 


Ap — AuJs{Pip — Pis) + AuJi{Pip — Pu) — 27P, (6) 

with Pim = dp/dujm = oiie over the group veloc¬ 

ity at each frequency, m = p^s^i. It can be seen that 
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the differences in ^jvg between pump, signal, and idler 
are an important factor in determining the JSA and its 
degree of correlation, or its factorability. Fig. shows 
three JSAs calculated approximately in this fashion. In 
Fig. a), the group velocities are chosen arbitrarily with 
Pip > Pii > Pisi resulting in a highly correlated JSA. 
In Fig. [^b), the idler is group velocity matched to the 
pump, resulting in an uncorrelated JSA, as the main peak 
now has its axes horizontal and vertical, although the 
ripples to either side, resulting from the sine function in 
G, remain correlated. This is the asymmetric scheme 
to generate a factorable JSA m- In Fig. B c), I3is and 
Pli are roughly equally spaced above and below /3ip, the 
symmetric scheme for a factorable JSA [16]. Here, the 
bandwidth of the pump must be exactly tuned to make 
the main peak circular in shape and uncorrelated. The 
conditions on the group velocities are generally met by 
using birefringent phase-matching, with the pump polar¬ 
ized on the slow birefringent axis and the photons on the 
fast axis. 

For a given JSA, the degree of correlation can be cal¬ 
culated using the singlular value decomposition function 
in Matlab. This provides a Schmidt decomposition of the 
JSA: 


3SA = '£^jfji^s)9M), ( 7 ) 

3 

where the fj{uJs) and gj{uJi) are a set of orthogonal spec¬ 
tral modes for signal and idler, and the Xj are real am¬ 
plitude coefficients. We define the purity as 

p = (8) 

3 

which is an upper limit on the quantum interference 
visibility possible between two photons from separate 
pairs m- For the JSAs in Fig. the purities are found 
to be 23% for the correlated JSA in (a), 95% for the 
asymmetric uncorrelated case in (b), and 83% for the 
symmetric case in (c). In (c), the purity is mainly limited 
by the correlation in the sine function ripples to either 
side of the main peak. These can be removed by filter¬ 
ing, resulting in high purities with relatively little cost to 
transmission efficiency m- 

The nonlinear correction —2^P to A/d, which has been 
ignored in the calculations above, depends on the pump 
power. For a pulsed pump, P is a function of position and 
time, and cannot be simply expressed in the frequency 
domain- previously the peak power has been used. In¬ 
stead, we take SPM and XPM into account by working 
in the time domain, and developing expressions for the 
joint temporal amplitude (JTA), which is linked to the 
JSA by 2D fourier transform. 


EQUATIONS OF MOTION IN A NONLINEAR 
FIBRE 


To model photon pair production through FWM with 
SPM and XPM included exactly, we use the equations 
of motions with position and time for the pump, signal, 
and idler fields in a medium. The electric field as¬ 
sociated with the pump pulse can be split into positive 
and negative frequency components as: 

E = eX{x,y) ^ ^-^iPpoz-ujpot)^ ^ 

(9) 

Here, e is the polarization vector of the electric field. 
X{x^y) is the transverse mode shape, normalised such 
that JJ X{x^y)‘^dxdy = 1. Complex oscillatory terms 
have been separated out at a central frequency ujpo and 
wave vector PpQ. This leaves £’+ as a complex envelope 
function describing the pulse, dependent on time t and 
position 2 : along the propagation axis 2 :. E~ is the com¬ 
plex conjugate of £’+. 

Assuming the pulse envelope varies slowly (the length 
of the pulse contains many optical cycles, or equivalently 
the bandwidth of interest is small compared to the fre¬ 
quency co’po), and that the nonlinearity is a small pertur¬ 
bation to the linear evolution of the pulse, obeys the 
nonlinear Schrodinger equation [21]: 


dz dt 2 dP 


il,\E+\^E+. 


( 10 ) 


The Pip term corresponds to the group velocity of the 
pump pulse. In the following it is removed from the equa¬ 
tion as we consider all quantities in a moving reference 
frame. The P 2 p term gives rise to dispersion - higher 
order dispersion terms are neglected here, but can be in¬ 
cluded as higher-order time derivatives. The term on the 
right-hand side is the nonlinearity associated with SPM 
of the pump, with Xp an effective nonlinear coefficient: 


2cnmA ’ 

where rim is the refractive index at the frequency ujm and 
A the effective area of the transverse mode. 

Analagous equations of motion for the signal and idler 
fields can be written 


dE+ dE+ 02, d^E+ 

dz dt 2 dP 


0s(2\E+\^Et + E;^E-) 

( 12 ) 


dEt 

dz 


EPx 


dE+ 02i d^E+ 
dt 2 dP 


= i'n{2\E+\^E++E+^E;) 

(13) 

In the moving reference frame. Pis and Pu are taken to 
be group velocity terms relative to the pump (ie. Pis 
Pis — Pip and Pii Pii — Pip). The first term on the 
right hand side of these equations is XPM, as the pump 
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modifies the refractive index experienced by signal and 
idler. The second term relates to FWM, with the strong 
pump field creating a coupling between signal and idler 
fields. For simplicity, the central frequencies of the signal 
and idler have been chosen to be a point of exact phase¬ 
matching with the central frequency of the pump, so that: 


2cjpo ^sO ^io — 0 n 4i 

WpO - PsO - Ao = 0. ^ ^ 

Note that since the pump field is many orders of magni¬ 
tude brighter than the signal and idler, which on aver¬ 
age will contain less than one photon, terms representing 
SPM of signal and idler, and XPM from signal or idler to 
other fields, are ignored. Similarly depletion of the pump 
due to FWM is neglected. 

The coefficient will generally be three times 

smaller in nonlinear effects coupling fields of orthogo¬ 
nal polarization compared to fields which are all co¬ 
polarized [21]. In the birefringent phase-matching 
schemes considered, signal and idler are orthogonally po¬ 
larized to the pump, so this is incorporated by reducing 
7s and 7^ by a factor of three, while keeping 7^ as above. 

Although the pump laser can continue to be treated 
classically, the signal and idler fields should be quantised. 
We use the following quantisation, similar to [22l [23] : 

where the a^j are annihilation operators for a photon at 
position z with frequency cj, with Hermitian conjugate 
creation operators dj,. An identical expression applies 
to the idler, with the integral taken to be over a dif¬ 
ferent range of frequencies, so that the two remain dis¬ 
tinct. The frequency modes are continuous, so the cre¬ 
ation and annihilation operators have a commutation re¬ 
lation [de^, dj^,] = — and the number density oper¬ 

ator dj,de^ should be integrated over a frequency interval 
in order to refer to the actual number of photons within 
that interval. 

It is convenient to consider the signal and idler in terms 
of the creation and annihilation operators for a photon at 
a particular time and position, (z, t) and t), which 
are the fourier transforms (from frequency to time) of 
alj and If the frequencies of interest for signal and 
idler he in a narrow bandwidth about and the 
electric fields have a simple approximate relation to the 
new operators: 


Et = 


huJsO 
-( 

2eocns 


E+ = 


huio 

-( 

2eocni 


( 16 ) 


is independent of signal and idler, unless we neglect the 
dispersion terms ^ 2 - Fortunately, as above, the purity 
largely depends on the different group velocities for the 
pump, signal, and idler, and the approximate solutions 
ignoring dispersion are still instructive. To include dis¬ 
persion properly, numerical methods can be used, as de¬ 
scribed later. 


APPROXIMATE SOLUTIONS 

Once the Pi term is removed from equation [Tq| using a 
moving reference frame, and the P 2 term is ignored, the 
pump pulse will retain its temporal shape as it propagates 
down the fibre, only accumulating a nonlinear phase due 
to SPM: 


E+{z) = E+{Q)e^<^- (17) 

where 

0p=7p^|^t(O)l' (18) 

Below, £’+ is taken to mean the pump amplitude as a 
function of time at 2: = 0. The signal and idler equations 
still include a group velocity term, in addition to XPM 
and FWM terms. In the interaction picture of quan¬ 
tum mechanics, the group velocity and XPM parts of 
the evolution, which affect signal and idler individually, 
are applied to the operators, while the FWM interaction 
between signal and idler is applied to the wavefunction 
IVi)- We first write the solutions ps^i to the group velocity 
and XPM terms, ignoring FWM: 

4>s{z,t) = 4>s{0,t-(3isz)E‘^‘, (19) 

where the nonlinear phase aquired due to XPM is 



with an analogous solution for pi. Note that because the 
signal (or idler) experiences group-velocity walk-off from 
the pump it accumulates phase from the pump at a range 
of different times, hence the integral. If the signal were 
group velocity matched to the pump. Pis = 0, the phase 
would become Os = 2^sz\E^{t)p. 

The evolution of the wavefunction according to FWM 
is now given by: 

^|V’)=iff|^), (21) 

with 


Like EE^ the ps,i have had the quickly varying oscilla¬ 
tions with 2: removed. 

These equations of motions do not have convenient so¬ 
lutions, even the classical equation for the pump, which 


H = I dtE+^<i>lcp\e^^'^^ + ( 22 ) 

Since the pair rate per pulse will generally be small, to 
avoid multi-pair emission, we take the interaction to first 
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order, beginning with the signal and idler modes in the 
vacuum state: 

l-tp) = \vac) + i^ 7 s 7 j jj dtdzE^‘^4>l4>\e^'^^^ \vac ), 

( 23 ) 

with L the fibre length. To extract the JTA from this 
wavefunction, we take the overlap between \^Ij) and a 
signal-idler pair at times 

JTA{ts,ti) = {vac\(t)s{L,ts)(t)i{L,ti)\i;), ( 24 ) 

Substituting in equation and simplifying we have: 


because there is no walk-off between pump and idler. 

Figure shows the joint temporal amplitude for three 
different lengths of fibre, with a Gaussian pump shape 
(X e ^ . It can be seen that for short lengths, the 

photons are highly correlated in time, with hard edges 
to the JTA due to the condition 0 < < A. As the 

length is increased, the signal walk-off smears out the 
JTA, making it closer to rectangular and less correlated. 
The fibre length is not important in itself so much as 
the ratio between the fibre length and the length of the 
pump pulse r, as changing both by a constant factor is 




0 


if 0 < Zc < T 
otherwise 


( 25 ) 


with Zc the point in the fibre at which the pair was cre¬ 
ated. This is defined for a particular tg^ti because the 
signal and idler must be created at the same time, G, 
and the extent to which they have walked off from each 
other identifies the length they have propagated through 
after creation. Similarly tc is defined by the differing 
arrival times of signal and idler: 


Zc = L - 


i^ls - Pli 


Pis - Pli 


( 26 ) 


0 is the nonlinear phase containing the contributions 
from SPM of the pump up until the pair-creation, and 
XPM from the point of pair-creation until the end of the 
fibre: 


e = \E+fdt+^^ 

( 27 ) 

The total probability of generating a pair, or the gen¬ 
eration rate per laser pulse, is given by 


R = 11 \JTA\^dtsdti = J ( 28 ) 


Clearly this solution becomes unphysical if Pig = Pu^ in 
which case the integrals implicit in equation!^ need to be 
dealt with differently. Also when dispersion is included 
it will affect R by causing the pump pulse to broaden or 
compress in time, respectively decreasing or increasing 
the generation rate. In the following R is by numerical 
integration of a JTA over tg and 


Asymmetric scheme 



When the idler and pump are group velocity matched, 
Pii = 0, tc = ti. Also the last term in 0, representing 
XPM from pump to idler, becomes 

( 29 ) 


FIG. 2. Joint temporal amplitudes when the idler and pump 
are group velocity matched. For increasing fibre length com¬ 
pared to the length of the pump pulse r, the group velocity 
walk-off of the signal smears out the JTA. (a) fSisLjr = 2 (b) 
/3isL/r = 5 (c) /3isL/r = 10. 


2yi{L-z,)\E+{U)\^ 














6 





FIG. 4. Purity plotted against generation rate for four choices 
of fibre parameters. Red: I3isL/t = 40. Green: ^igLlr — 
20. Blue: ^xsLjr — 10. Purple: PisL/r — 5. The black 
dashed line shows the purity varying with generation rate 
after numerical optimization oi ^\sL j r . It can be seen a large 
value of PisL/t is desirable at low rates, but for higher rates 
it will cause a more rapid fall-off in purity. 



FIG. 3. JSA for different generation rates R, showing the 
effect of SPM and XPM. The initial effect is to broaden the 
idler, eventually leading to a splitting and distortion of the 
JSA which reduces the purity P. (a) R approaching zero, 
P = 89% (b) R = 0.1, P = 83% (c) R = 0.2, P = 78%. 

simply a rescaling of the JTA. It will appear uncorrelated 
if PisL/t ^ 1, though in practice this may be limited by 
dispersive effects which have been ignored so far, since 
they will become more significant for longer lengths and 
shorter pulses. 

The JTAs are shown as absolute values and so are not 
affected by the nonlinear phase. The effects can be seen 
in the JSA obtained by taking the fourier transform of 
the JTA - in figure the JSA is shown for increasing 
probability of pair-creation in a pulse, in each case 
with PiqL/t = 10. As R increases, the idler is broadened 
significantly by phase modulation and begins to distort 


in profile. This causes the purity to decrease, from 89% 
at = 0, to 83% at = 0.1, to 78% at = 0.2 (0.2 pairs 
per pulse may be unrealistically high for an experiment, 
but allows the distortion of the JSA to be seen clearly). 

In Figure]^ the purity is plotted against the pair gen¬ 
eration rate R for different choices of fibre parameters. It 
can be seen that for larger values of pigL/r^ the purity 
will be high at very low i7, but will decrease rapidly as R 
increases, whereas a smaller value oi Pi qLIt will decrease 
more slowly and may be optimal for a given R. Note that, 
even after optimising the purity with the fibre parame¬ 
ters, the purity will be more detrimental to quantum in¬ 
terference quality than multi-pair emission over the range 
shown, 0 < < 0.1. [The probability of multi-pair emis¬ 

sion is estimated as which when compared to the 
rate of single-pair emission R can reduce the interference 
visibility by at most R]. This is potentially significant 
if the end goal is to build a deterministic photon source 
by multiplexing together many of these sources [9] , then 
to achieve high-quality interference without filtering for 
quantum communication or computing applications. It is 
usually assumed that the end quality will be high so long 
as multi-pair emission is kept low from the individual 
sources, but this shows that, at least for this asymmetric 
scheme using FWM, the effects of phase modulation are 
likely to be the limiting factor on the generation rate. 

Inspection of equation does suggest a solution to 
this problem. The nonlinear phase factor e*® becomes a 
factorable function of tg and U over the extent of the JTA 
so long as the pump field is a square function in time. 
The phase is only correlated because \E^{tc)\‘^ is varying 
across the JTA. However, the effects of group-velocity 
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dispersion acting on a short, square pulse over a large 
length may be unpleasant. Using realistic dispersion pa¬ 
rameters based on the birefringent microstructured fibre 
in [24l|25] with a length 50cm, and calculating the JSA 
from equation without phase modulation effects, the 
maximum value of P using a square pulse of optimal du¬ 
ration is found to be 80%. Prechirping the pulse to com¬ 
pensate the dispersion, so that it is square at the mid¬ 
point of the fibre, yields a slight improvement to 81%. 
So this is unlikely to be helpful unless the dispersion is 
particularly small. 

Another solution would be to have 7^ = 7i, although 
this is not possible using the birefringent phase-matching 
considered, because of the reduction in the effective non¬ 
linearity by a factor of 3 when the fields are orthogonally 
polarized. However if the nonlinear phase could be made 
factorable, a high purity could again be achieved with a 
large value of pigL/r. 

Symmetric scheme 

We now consider the symmetric scheme for avoiding 
correlations, with Pis = —Pu = Pi. This implies that 
Zc = L — and tc = Figure ^ a) shows the 

JTA in this case, where the temporal wiom r of a gaus- 
sian pump has been optimised to minimise correlation. 
Figure |^b) shows the corresponding JSA without the 
effects of phase modulation, with R approaching zero, 
while Figure |^c) shows the broadening and distortion 
from phase modulation when R = 0.2. Here, the broad¬ 
ening introduces spectral correlation, but it can be partly 
compensated by beginning with a longer pump pulse (in¬ 
creasing r). Figure shows the purity plotted against 
generation rate, both for a fixed value of PiL/r^ and 
with r reoptimised as R is increased. 

The predicted purities of around 80% here are some¬ 
what low, even when R is kept small. A realistic ex¬ 
perimental strategy may be to increase the purity above 
this by applying some spectral filtering to one photon of 
the pair, with an overall transmission T. If the filtered 
photons are used as the heralds, the heralded photons 
will not experience any loss, just a reduction in effective 
generation rate to RT. Figure shows the purity af¬ 
ter one of the photons has been filtered with a top-hat 
transmission window of variable width, as a function of 
RT. Five different values for the original generation R 
are shown, from 0.2 to 0.1. As expected the purity tends 
to 1 as the filtering becomes very drastic and only leaves 
one possible spectrum for the heralded photon, although 
this restricts the source to low effective generation rates. 
It can also be seen that the purity tends to 1 faster if the 
initial (unfiltered) generation rate R is larger, in spite of 
the detrimental effects of phase modulation. 

Since a larger R appears to be beneficial for the purity 
after filtering, as a function of i?T, this suggests there 



FIG. 5. (a) Joint temporal amplitude when the signal and 
idler are equally spaced in Pi about the pump. The pump 
duration is optimised to avoid correlation, but the hard edges 
to the JTA caused by the sudden beginning and end of the 
nonlinearity mean that some correlation is inevitable, (b) 
Corresponding JSA with low generation rate, (c) At higher 
generation rate, R = 0.2, showing the effects of SPM and 
XPM. The JSA is broadened in one direction. 


will be a trade-off between achieving higher purity and 
keeping multi-pair emission low, which occurs with prob¬ 
ability approximately R^ m- If photon-number resolv¬ 
ing detectors become available, they could be used on the 
heralds to detect and filter out multi-pairs. Otherwise it 
may not be possible to simulataneously achieve a high 
purity, high effective generation rate, and low multi-pair 
emission using this scheme. In future work, it would be 
useful to consider similar schemes where two pump fre- 













FIG. 6. Purity against generation rate for the symmetric 
scheme. Blue: r is kept constant. Green: r is increased with 
R to reoptimise. Again, the phase modulation reduces P as 
R is increased, though here P is lower to start with than in 
the asymmetric scheme because there is more correlation in 
the sine ripples of the JSA. 
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FIG. 7. Purity as a function of the effective generation rate, 
RT, when hltering is applied to one of the photons (the her¬ 
ald) resulting in transmission T. Unhltered generation rates 
R are shown from 0.02 to 0.1. 


quencies are used, with mismatched group velocities m- 
This can in principle remove the sine ripples from the 
JSA, and hence most of the correlation, although phase 
modulation may still have an effect. 


NUMERICAL MODEL 


To include the effects of group-velocity dispersion accu¬ 
rately, it is necessary to go to a numerical model involving 
finite-steps along the fibre length. A common method for 
modelling the propagation of a laser pulse through a non¬ 


linear and dispersive medium is a split-step fourier (SSF) 
simulation [21]. Here, the length is divided into small 
steps Az, and for each step, the nonlinearity and the dis¬ 
persion are applied separately. For instance, the effect of 
propagating through the nonlinearity of Az could be ap¬ 
plied first, in the time domain where this is a simple cal¬ 
culation, then the pulse could be fourier transformed to 
the frequency domain, where the effect of the dispersion 
can easily be applied using E~^{z, Oj) = E+{0, 
followed by inverse fourier transform back to the time do¬ 
main. Since the effects of the nonlinearity and dispersion 
are generally non-commuting, this is only approximate, 
but is accurate for small Az. In fact, it is better to apply 
half a step of dispersion, then a full step of nonlinearity, 
followed by the other half step of dispersion, as then the 
size of the errors due to the approximation vary with Az^ 
rather than Az‘^ m- 

To model the pair-production process along similar 
lines, we use a SSF simulation for the propagation of 
the pump pulse, and incorporate spontaneous FWM into 
the nonlinear part of each step. The steps in position 
are kept small compared to the resolution in time, so 
that Pis^iAz < At. This means the signal and idler are 
intially in the same time-bin as each other, and as the 
component of the pump which created them, which sim¬ 
plifies the calculation for each step. To propagate the 
state of pairs created in previous steps, the half-step of 
dispersion is applied to the JSA, which is then converted 
to a JTA by 2D fourier transform so that XPM can be 
applied, before it is transformed back to a JSA for the 
next half-step of dispersion. The state of the new pairs 
is coherently added for each step. 

In the previous sections, the purity was determined by 
two parameters: the group velocities relative to the fibre 
length and pulse duration, PiL/r, and, when nonlinear 
phase modulation was taken into account, the total prob¬ 
ability of pair-generation, R. Here, the group-velocity 
dispersions at each wavelength introduce additional rel¬ 
evant parameters: /32L/r^ for pump, signal, and idler. If 
the fibre length and pulse duration are increased in pro¬ 
portion the effect of group-velocity dispersion decreases. 

We again consider dispersion parameters taken from 
the birefringent microstructured fibre used in [22 1 ^ . 
This fibre makes use of the asymmetric scheme to avoid 
correlations, with the pump pulse at 72Qnm group- 
velocity matched to the idler at 864nm. The sig¬ 
nal, phase-matched at 626nm, experiences walk-off with 
Pis = 1.14 X Pump, signal, and idler ex¬ 

perience group-velocity dispersion with P 2 p = 2.1 x 
P2s = 3.6 X and P2i = —1.3 x 

. Figure shows the JSA produced from a 
50cm fibre with a 2nm initial pump bandwidth, corre¬ 
sponding to r ^ 230/s, and with the generation prob¬ 
ability R = 0.05. The dispersion introduces some cur¬ 
vature to the JSA, so that it will become correlated for 
large bandwidth pulses. 
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FIG. 8. JSA from the numerical model including dispersion, 
for realistic fibre parameters with a length of 50cm and a 
pump bandwidth 2nm. The dispersion causes some curvature 
of the JSA, which may introduce correlation. 


In figure the purity is plotted against the gener¬ 
ation probability using the numerical model, for differ¬ 
ent amounts of dispersion: no dispersion, dispersion us¬ 
ing the realistic parameters, and double strength disper¬ 
sion. Physically, doubling the strength of the dispersion 
while keeping the other parameters constant could be 
achieved by halving the length of the fibre, halving the 
laser pulse duration, and adjusting the laser power to 
keep R constant. In figure [^a), I3isL/t is approximately 
12. It can be seen that for a very low generation rate, 
where phase modulation is negligible, the dispersion only 
causes a slight reduction in purity. The zero R intercept 
is close to our previous modelling results which ignore 
non-linear phase modulation [15]. However, as R is in¬ 
creased, dispersion has a larger effect, suggesting that 
dispersion and phase modulation are combining to cre¬ 
ate more of a reduction in purity than either would alone. 
Surprisingly, for higher values of i?, this trend reverses, 
and the case with some dispersion actually does better 
than no dispersion. This can be seen more prominently 
in figure [^b), where the fibre length was doubled to Im, 
so that I3isL/t ^ 2A. It seems that for some choices of 
parameters, the nonlinear phase modulation and the dis¬ 
persion begin to compensate one another, although here 
the effect is too small to change the trends seen before, 
with the purity still decreasing as the generation rate is 
increased, and with a larger value of PigL/r creating a 
better purity at low rate, but a worse purity at high rates. 

Finally, we consider a different set of fibre parame¬ 
ters, corresponding to a birefringent microstructured fi¬ 
bre pumped at 1064nm, in its anomalous dispersion re¬ 
gion, generating phase matched photons at 810nm and 
1550nm m- The pump is polarized on the fast axis 
of the fibre while signal and idler are polarized on the 


slow axis, such that the signal is now group velocity 
matched to the pump, while the idler walks off with 
Pii = 1.2 X The dispersion parameters are 

/32p = —8.7 X ^23 = 1-0 x and 

hi = —6.4 X Figure [lQ|^a) shows the purity 

plotted against generation rate for a 50cm length of this 
fibre with an initial pump bandwidth of 2nm, resulting 
in hsL/r 12, and figure [l^b) shows the case when 
the length is increased to Im, so hsL/r ~ 24. It can 
be seen that at low i?, the dispersion improves the purity 
slightly compared to the case without dispersion but with 
phase modulation. However, for larger generation rates 
the combination of dispersion and phase modulation has 
a significant deleterious effect, as can be seen clearly in 
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FIG. 9. Purity against generation rate R for different 
strengths of dispersion, with an initial pump bandwidth of 
Inm. R — 0 corresponds to the case with no non-linear phase 
modulation, (a) L — 0.5?n,(b) L = Im. In each case, black 
circles: no dispersion; blue squares: realistic dispersion; red 
triangles: double strength dispersion. Surprisingly, at higher 
rates the dispersive case sometimes does better than the case 
with no dispersion. 
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FIG. 10. Purity against generation rate with anomalous dis¬ 
persion, for an initial pump bandwidth of 2nm. (a) L = 
0.5m,(b) L = Im. In each case, black circles: no dispersion; 
blue squares: realistic dispersion. Here, a very slight benefit 
can be seen at low rate for the dispersive case compared to 
the case without dispersion. 


figure [T^b). 


CONCLUSION 

We have seen that for schemes seeking to minimise the 
correlation between photon pairs generated by four-wave 
mixing, the effects of self-phase modulation and cross¬ 
phase modulation may be a limiting factor on the pho¬ 
tons’ purity which is not usually considered. For the 
asymmetric scheme, where one of the generated photons 
is group velocity matched to the pump pulse, it would 
otherwise be expected that, with a long interaction length 
and a wide pump bandwidth, very high purity could be 
achieved. However, when these additional nonlinear ef¬ 
fects are included, the purity is degraded as the gener¬ 


ation rate is increased, which may limit sources to low 
rates when a particular purity or quantum interference 
visibility is required. This can be seen both in an ana¬ 
lytical model in the time domain, where group-velocity 
dispersion is neglected, and in a numerical model which 
includes both nonlinear effects and dispersion. 

The symmetric scheme to generate pure photons is 
also considered, with the signal and idler group indices 
equally spaced above and below the pump group index. 
Again the nonlinear effects significantly degrade the pu¬ 
rity as the generation rate is increased, although here the 
main source of impurity is the correlation in the sinc- 
ripples of the phase-matching function. These ripples 
can be eliminated with narrow filtering, but the purity 
only tends to unity as the transmission through the fil¬ 
ter becomes small. In future work it would be interest¬ 
ing to consider the case with two pump fields at different 
wavelengths, where the sine ripples can in theory be elim¬ 
inated without filtering, but it seems likely that nonlinear 
effects will also be detrimental there. 

The numerical model demonstrates that the impu¬ 
rity from nonlinear effects and from dispersion do not 
combine trivially, sometimes leaving a lower purity than 
would be expected when the effects are taken individu¬ 
ally, but in some regimes slightly higher. It is possible 
that for particular choices of pump pulse power, dura¬ 
tion, and shape, the effects of dispersion and nonlinear¬ 
ity could be made to cancel in a soliton-like manner, and 
leave a high purity, although it is not expected that hav¬ 
ing the pump pulse alone propagating as a soliton would 
achieve this. 

We conclude that the discrepancy between previous 
modelling results and measured visibilities in heralded 
photon interference experiments nails] can be largely 
explained by including non-linear phase modulation. We 
also note that this effect limits the fidelity of cluster states 
generated by fusing entangled states [24l [25] and thus 
could limit the scalability of cluster state quantum com¬ 
putation based on four wave mixing 
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